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Abstract 

Dynamic simulations and experimental validation tests 
were performed on a two-stage, two-speed gearbox as part of 
the drive system research activities of the NASA Fundamental 
Aeronautics Subsonics Rotary Wing Project. The gearbox was 
driven by two electromagnetic motors and had two 
electromagnetic, multi-disk clutches to control output speed. A 
dynamic model of the system was created which included a 
DC electric motor with proportional-integral-derivative (PID) 
speed control, a two-speed gearbox with dual electro- 
magnetically actuated clutches, and an eddy current 
dynamometer. A six degree-of-freedom model of the gearbox 
accounted for the system torsional dynamics and included 
gear, clutch, shaft, and load inertias as well as shaft 
flexibilities and a dry clutch stick-slip friction model. 
Experimental validation tests were performed on the gearbox 
in the NASA Glenn gear noise test facility. Gearbox output 
speed and torque as well as drive motor speed and current 
were compared to those from the analytical predictions. The 
experiments correlate very well with the predictions, thus 
validating the dynamic simulation methodologies. 
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4 

c 
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fni[i =1,2] 

.//;/, max 

F c , F tl, F Tm 
ho, hd 

Ic 


EMA pole face area 
damping matrix 
set-point speed error 
EMA clutch actuation force 
maximum actuation force 
generalized forces 
coil static and dynamic air gap 
EMA coil current 


Umax, Usal 

maximum and saturation coil current 

U: 

motor armature current 

i 

i th clutch or i th gear 

J cdi, J cfi 1 t 1 -2 | 

clutch disk and clutch flange inertias 

Jl, Jm 

dynamometer and motor inertias 

J gi [i 1,6] 

gear rotational inertias 

kp, ki, kd 

motor PID speed control gains 

k,nc\, k c \ c l, k s l , k c 2L 

shaft torsional stiffness's 

kmb, knit 

motor back EMF and torque constant 

k s 

clutch separator spring constant 

K 

stiffness matrix 

km 

motor circuit inductance 

M 

mass matrix 

n n, W 34 , w 56 

gear ratios 

N d 

number of clutch disks 

N w 

number EMA coil windings 

q 

DOF vector 

q/7, q hi, q iv 

constrained DOF vectors 

Q 

total generalized force vector 

R/n, Rc 

motor and EMA coil electrical resistance 

Rd, Ro 

clutch disk inner and outer radii 

t 

time 

T 

system kinetic energy 

Tl, T m 

dynamometer and motor torques 

T ci [i= 1,2] 

EMA clutch torques 

T c i,n, T ci ju, T ci i V 

clutch torques (kinematically locked) 

Tkmax U = 1,2] 

maximum static clutch torques 

Tci.slip[i= 1,2] 

slipping clutch torques 

V 

system strain energy 

V ci [i = 1,2] 

EMA clutch coil voltage 

V 

v cmax 

maximum clutch coil voltage 

v d 

dynamometer voltage 

v m 

motor armature voltage 

Zip =1,6] 

number of gear teeth 
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Stribeck exponent (empirical) 
virtual work 

slip-stick threshold tolerance 
clutch disk and flange rotation 
gear rotations 
load and motor rotation 
free space permeability (air) 
kinetic and static friction coefficients 
Dry friction coefficient (Stribeck) 
viscous damping parameter 
clutch disk and flange speed 
load and motor speed 
motor set-point speed 
slip speed 

Stribeck velocity (empirical) 


Introduction 

As part of the Subsonic Rotary Wing Project of the NASA 
Fundamental Aeronautics Program, research is being 
performed for large civil tiltrotors to replace regional airliners 
over medium ranges to alleviate next-generation air traffic. 
The tiltrotor is designed to carry 90 passengers for 1000 nm 
with performance of 300 knots at 28,000-ft altitude. Aircraft 
sizing-design studies have identified that a two-speed rotor 
configuration is required (Refs. 1 and 2). Rotor speed at cruise 
is required to be reduced, possibly down to 50 percent of 
hover speed, to keep rotor blade advancing tip speeds 
reasonable. To support this configuration, research in 
developing two-speed gearboxes as well as dynamic modeling 
of multi-speed rotorcraft drive systems is being led at the 
Glenn Research Center with help from university counterparts. 

As a brief background, one early study of a multi-speed 
rotorcraft drive system consisted of a high-speed traction drive 
variator and a planetary differential gear unit (Ref. 3). A unit 
was tested and performed as designed. However, the 
mechanical efficiency was lower than that of a conventional 
drive system. Also, the use of traction drives for a main drive 
path is not an accepted practice for U.S. rotorcraft 
manufacturers. Previous studies have been sponsored by 
NASA to investigate multi-speed drive systems (Refs. 4 and 
5). These studies looked at earlier tiltrotor applications where 
15 percent speed reduction for cruise was required. Dual-path 
configurations and compound planetary configurations looked 
promising. However, major concerns that were found included 
shifting and drive system weight. An electro-mechanical, 
infinitely-variable transmission, comprising a pair of planetary 
trains interconnected with two electric machines and clutches, 
has been proposed (Ref. 6). Again, the mechanical efficiency 
suffered as well as the added complexity of two planetary 
systems. Also, a unique concept called a pericyclic 
continuously variable-speed transmission is being investigated 
(Ref. 7). It uses pericyclic kinematics to achieve single-speed 
reduction ratios between 1.05:1 to 50:1 with variable-speed 


capability in one configuration. This concept, however, has yet 
to be tested. In summary, the results from the previous studies 
indicates that the incorporation of multi-speed concepts in 
rotorcraft application is not a trivial process and adds 
complexity and weight. Thus, further research and 
development in this area is required, such as described in 
(Ref. 8). in addition, dynamic modeling of multi-speed drive 
systems can help define the characteristics and limitations of 
the drive systems and the shifting process. 

Since very few multi-speed rotorcraft studies have been 
published, very little work on dynamic modeling for multi- 
speed rotorcraft drive systems has been done. There has been, 
however, numerous publications on dynamic modeling of 
automotive drive systems. Much of the work addresses the 
automatic transmission shifting control with emphasis on 
maximizing fuel economy and maintaining performance 
(Refs. 9 to 15). For these, the typical system being modeled is 
the gasoline internal combustion engine, torque convertor, 
multi-speed gearbox, tires and vehicle dynamics, as well as 
throttle and shifting controls. Dynamic speed, gear ratio, and 
torque are some of the parameters predicted. This strategy has 
been applied to hybrid power trains (Refs. 16 and 17), 
continuously variable transmissions (Refs. 1 8 to 22), and even 
heavy trucks and tanks for the military (Refs. 23 to 27). The 
modeling in the automotive field has had much success in 
improving drive train development. The approach and some 
modeling tools can be extended to the rotorcraft arena. 
However, specifics of integration of the gas turbine engine, 
multi-speed drive systems, clutches, and the dynamics of the 
rotor and rotorcraft vehicle have yet to been addressed for 
rotorcraft application. 

A NASA Research Announcement (NRA) contract was 
awarded to the team of Penn State University, the University of 
Tennessee, and the University of Michigan, to develop a 
dynamic model of a multi-speed rotorcraft drive system. The 
model integrates the dynamics of the gas turbine engine, multi- 
speed transmission system, and rotor system. System models for 
conventional helicopters and tiltrotors are being developed. 

The objective of the current study is to validate a gearbox 
dynamic model developed under the NRA award. As a first 
step in the model validation process, a dynamic model of a 
two-speed gearbox developed at NASA for lunar space 
application was derived. This was an existing two-speed 
gearbox available for test. Experimental tests were performed 
on the gearbox as installed on a dynamometer. Speeds and 
torques were measured and compared to predictions during 
various shifting profiles. 


Chariot Gearbox Description 

The Glenn Research Center designed and fabricated a gearbox 
for use on the NASA Chariot vehicle. The Chariot is the next- 
generation rover vehicle for lunar exploration (Fig. 1) and was 
designed by the NASA Johnson Space Center (JSC) for truck- 
type utility operations on the moon. This includes the 
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\ ' Steering gearbox 
' Differential gearbox 

Figure 1 . — NASA Chariot lunar rover vehicle. 


ability to carry astronauts and substantial amounts of 
equipment, rocks, or soil. It can be fitted with a soil-moving 
blade and the flat-top deck of the vehicle lends itself to custom 
configurations for different missions. It can also be outfitted 
with a pressurized crew cabin as shown in Figure 1. 

The Chariot has six identical and independently-controlled 
motor-driven wheel pods, three on each side of the vehicle. 
Each pod has a pair of driven wheels for a total of 12 wheels. 
Multiple wheels and pods creates redundancy for reliability and 
reduces the tire surface contact pressure. All six pods and wheel 
pairs can be steered independently and to any angle (a hill 360° 
capability), making the vehicle highly maneuverable. Each pod 
has its own combination passive and active suspension. The 
active part allows for adjusting the chassis height, including 
lowering it to the ground, which allows for easy egress of the 
astronauts and easier loading and unloading of cargo. 

Each pod is driven by a motor/gearbox assembly. The 
output passes through a steering gearbox, then to a differential 
gearbox, and then to each wheel (Fig. 1). The motor/gearbox 
assembly is the subject of the current study and will be further 
called the Chariot gearbox for simplicity. The original Chariot 
gearbox was designed by JSC and is shown in Figure 2. Due 
to a tight schedule, the JSC Chariot gearbox was designed and 
fabricated using readily available off-the-shelf commercial 
gears, bearings, shafts, electric clutches, and electric motors. 
Two brushless electric servo motors drive the unit, both to 
generate adequate power and for motor redundancy. The 
motors are rated for 300 V DC and each one draws up to 8 A 
to produce 2.8 hp at peak efficiency. The original gearbox 
design utilizes a total of 12 gears, 17 bearings, and 7 different 
shafts. It has two selectable speed ratios. The low speed 
provides a 16:1 gear reduction for a vehicle speed of about 
3 mph while the high speed provides a 4: 1 gear reduction for a 
vehicle top speed of 12 mph. The speed selection is 
accomplished with two electric clutches. 

The NASA Glenn Research Center (GRC) used their 
rotorcraft drive system experience to redesign the Chariot 



Figure 2. — NASA Johnson Space Center Chariot gearbox 
design. 


gearbox. The gearbox was redesigned as a direct replacement 
to the JSC design, but used gear design practices of the current 
rotorcraft industry to reduce the number of components as 
well as weight and size. To allow direct replacement in the 
Chariot vehicle, the same servo motors were used in the 
redesign as well as the same gearbox mounting pattern and 
output shaft connection. Custom-made high-strength gears 
allowed the use of lighter and fewer gears (7 instead of 12), 
and the gear train consists of just two parallel axes, resulting 
in fewer shafts (3 instead of 7) and fewer bearings (13 instead 
of 17). A detailed description of the redesign is given in 
(Ref. 28). 

A schematic of the GRC Chariot gearbox is shown in 
Figure 3. The gearbox is a two-stage, two-speed gearbox. The 
outputs of the two brushless electric servo motors are attached 
to pinion gears (Gear 1 in Fig. 3). For the first reduction stage, 
the pinion gears drive a bull gear (Gear 2) producing a 3.94:1 
speed reduction. The bull gear (Gear 2) drives two concentric 
electromagnetically-actuated (EMA) multi-disk clutches. For 
low-speed operation, the top clutch (low-speed clutch) is 
engaged and the bottom clutch (high-speed clutch) is 
disengaged. The output of the low-speed clutch passes through 
two sets of spur gears (Gears 3, 4, 5, and 6). This produces a 
16.14:1 speed reduction. For high-speed operation, the bottom 
clutch (high-speed clutch) is engaged and the top clutch (low- 
speed clutch) is disengaged. The bull gear (Gear 2) output 
passes directly out the gearbox producing a 3.94:1 speed 
reduction. Also, Gears 3 to 6 freewheel under no load during 
high-speed engagement. Design parameters of the gears are 
given in Table I. The gears were made from AISI 9310 steel 
and were carburized, hardened, and ground to an AGMA gear 
tolerance Class 12. 
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TABLE I.— NASA GLENN RESEARCH CENTER CHARIOT GEARBOX GEAR DESIGN DATA 



Gear 1 

Gear 2 

Gear 3 

Gear 4 

Gear 5 

Gear 6 

No. of teeth 

16 

63 

20 

42 

21 

41 

Diametral pitch (teeth/in.) 

8 

8 

10 

10 

10 

10 

Circular pitch (in.) 

0.393 

0.393 

0.314 

0.314 

0.314 

0.314 

Whole depth (in.) 

0.294 

0.294 

0.235 

0.235 

0.235 

0.235 

Addendum (in.) 

0.125 

0.125 

0.100 

0.100 

0.100 

0.100 

Chordal tooth thickness (in.) 

0.193 

0.193 

0.156 

0.156 

0.156 

0.156 

Pressure angle (deg) 

25 

25 

25 

25 

25 

25 

Pitch diameter (in.) 

2.000 

7.875 

2.000 

4.200 

2.100 

4.100 

Outside diameter (in.) 

2.250 

8.125 

2.200 

4.400 

2.300 

4.300 

Root fillet (in.) 

0.050 

0.050 

0.050 

0.050 

0.050 

0.050 

Tip relief (in.) 

0.0005 

0.0005 

0.0005 

0.0005 

0.0005 

0.0005 



Low speed 
clutch 

3 


speed 

clutch 


Figure 3. — NASA Glenn Research Center Chariot gearbox 
design. 


Commercially available off-the-shelf clutches were used in 
the redesign due to schedule constraints. Electric multiple disc 
clutches were used. They provided a maximum torque of 
70 lb-ft and were powered by a 24- V DC source, similar to the 
JSC Chariot application. 

The gearbox housing was manufactured from a 
combination of 6061 and 7075 aluminum alloys. A stacked- 
component approach was used in the design of the housing. A 
total of five parts made up the complete housing: 1) upper 
cover, 2) transition case, 3) deep case, 4) short case, and 5) 
lower cover. This helped keep the size profile minimized 
because it allowed the top gear plane, with the two motor 
pinions and bull gear, to be shrouded closely by way of the top 
cover before transitioning to the majority smaller profile 
required for the main body which encompasses the rest of the 
gears and both electric clutches. The gearbox was grease 
lubricated, using a grease similar to that used in helicopter tail 
rotor shaft bearings. Grease -retention shrouds were designed 
around all gears to hold the applied grease close to the gear 
teeth. 



Figure 4. — Dual-clutch gearbox torsional dynamic model. 


Chariot Gearbox Modeling 

To predict the shift response of the dual-clutch gearbox 
system, a torsional dynamics model including gear, clutch, 
shaft and load inertias, as well as shaft torsional flexibilities, 
and dry clutch stick-slip friction effects was developed. A 
schematic of the model is shown in Figure 4 and the 
kinematics is summarized in Equations (1) to (3). 

^gl _l l/m — — ^12^gl 

‘tyi = <l>g3> <t>g4 = ~ n 34§g3 (1) 

4^6 =_w 56 < l , g5> §cf2=§g6 

Here, <f),„, is the motor rotation, § cd \, (jyg, are the low-speed 
clutch (clutch no. 1) disk and flange rotations, § cd2 , § c fi, are the 
high-speed clutch (clutch no. 2) disk and flange rotations, 
())„, [/ =1-6] are the gear rotations, and (fy is the rotation of the 
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driven load. Also, the gear ratios n 12 , n 34 and n 56 in terms of 
the number of gear teeth Z, [i =1-6] are 

«I 2 = Zi /Z 2 , «34 = Z 3 /Z 4 , and « 56 = Z 5 /Z 6 (2) 

The engagement and disengagement of the EMA clutches is 
controlled by axial forces, f, \ and f n2 , applied normal to the 
clutch disks via electromagnetic actuation forces. 

Due to the kinematics, the only condition under which both 
clutches can be locked simultaneously corresponds to the 
zero-speed or fully-stalled condition. When both clutches are 
applied simultaneously the clutches act as a braking system. 
Thus, during normal operation only one clutch is fully 
engaged and locked (with full applied engagement force 
f ni = fni'tmx ) while the other clutch remains open (with zero 
applied engagement force,/,, = 0). In the low-speed condition 
clutch no. 1 is locked and clutch no. 2 is open and the 
transmission output speed, II/. is 


8W= 2T m 8^> m - r cl 5(<l> c/ i - <|> crfl ) 

... - T c2 8(§ c f2 ~ §cd2 ) _ Tl§§l 


where T„, is the applied motor torque, T c] and T c2 are the clutch 
friction torques between the disk and flange, and 7} is the 
dynamometer load torque on the output shaft. 

The equations of motion are obtained using Euler-Lagrange 
equations 


d_ dT_ 
dt <9q 


dT_ dV_ 

5q 5q 


(7) 


resulting in 


Mq + Cq + Kq = F £ .(q,q,/„i,/„ 2 ) + F ri +F rAf (8) 


with generalized coordinate vector 


Low ratio: fl L = - — — — D. 

Z 2 Z 4 Z fi 


(3a) 


0 [^m §cd2 §cf2 ^2.] (9) 


and in high-speed condition clutch no. 1 is open and clutch no. 
2 is locked 


High ratio: D. L = 



(3b) 


where is the motor speed. 

Based on the unconstrained condition with both clutches 
slipping, the system kinetic energy can be expressed as 


T = 2 


" (dm F *7gl F H\2 J g2 


2 

2 


+ — 


J cd 2 


••• + ~(^c/l +^g3 + ^34 J g4)§cf\ +~^¥cd2 ( 4 ) 

• • • + - (J C J2 + Jg6 + -^y)^>c/l + 

A «56 A 


where, J cdi , J cfl [i = 1,2], J gi [i =1-6], and J L are the 
rotational inertias of the motor, clutch disks, clutch flanges, 
gears and load, respectively. Furthermore, the total system 
strain energy is 


^ ^ [^mcl (*l ) g2 fycdl)^ F k c \ c 2 (§cd2 fycdl)*' 

■■■ + k c 2L($L ~§cf2) 2 +^1(^5 _< l , g4) 2 ] 


(5) 


where k mch k cIc2 , k s/ , and k c2 L are the torsional stiffness values 
of the motor-clutch 1 shaft, clutch 1 -clutch 2 shaft, the low- 
speed shaft, and outputs shaft, respectively (see Fig 4). Finally 
the virtual work due to the system torques is 


and with mass and stiffness matrices, M and K, defined as 


M = 


2 (J m + J g\ + n \2 J g2 ) 
0 
0 
0 
0 
0 

0 

0 

0 

J cd2 

0 


0 0 

Jcd\ 0 

2 

0 J cf\ + dg3 +H34~Jg4 

0 0 

0 0 

0 0 

0 


0 

0 

0 

0 

J cf2 + Jg6 + 


7g5 

2 

"56 


0 


0 

0 

0 

0 

J /. 


(10a) 


K = 


kmcl n \l 2 

kmc\ n ll 

0 

0 

0 

0 

k,nc\ n 12 

k/itc\ 2- k c \c2 

0 

-k c lc2 

0 

0 

0 

0 

ksl n 34 2 

0 

- A -," 34 

K sl 

"56 

0 

0 

- k c \c2 

0 

^clc2 

0 

0 

0 

0 

k , "34 
K sl 

"56 

0 

k j- ^ sl 

k c2L + Y 

«56 

~ /c c2L 

0 

0 

0 

0 

~k c 2L 

kc2L 


(10b) 


together with an assumed proportional damping matrix 

C=^K (10c) 

based on a structural viscous damping factor, q, = 0.00001 
(constant). Finally, the generalized force vectors obtained 
from the virtual work are 
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8fV = Q T 8q , Q = F Tm +F C +F TL (11a) 


with 


(N 


0 


0 ~ 

0 


Tcl 


0 

0 

, F c = 

-Tci 

, and F tl = 

0 

0 


Td 

0 

0 


-Tci 


0 

0 


0 


1 

1 


(lib) 


where F Tm , F TL , and F r are generalized forces due to T m , T,, 
and clutch friction torques T cl and T c2 . 

The DC electric motor torque, T m , is computed based on the 
armature current, /„„ and applied motor control voltage, V m , as 


dl m _ j - 

, ' ^m 1 m ' v m 

dt 

T —Jr J 

1 iTi rv mt- L m. 


(12a) 


where L m , R m , k m b, k mt are the motor circuit inductance, 
resistance, back EMF constant, and torque constant, 
respectively. Furthermore, V m is determined using a 
Proportional Integral Derivative (PID) speed controller 

V m = k p e + k d e + kjje dt (12b) 


which regulates Q„, to the set point speed fl sel based on the set- 
point speed error e = Cl set - Cl m and the control gains k p , k d , and 
k,. 

As described in the experimental test setup (see next 
section), the output shaft of the gearbox is connected to an 
eddy current dynamometer. In the test mns, the dynamometer 
is set to operate in so-called constant torque mode. Thus, in 
the analysis model, the resistive load torque, 7), is modeled as 

Tl = T L ,set tanh(F rf Q i ) (13) 


where, N d is the number of friction disks in each clutch pack, 
and Rj and R„ are the clutch disk inner and outer radii, 
respectively. Furthermore, p s * is computed using the so-called 
Stribeck dry friction model (Ref. 29). 

V-stb =fbt +(Ps -P/Oe“ ( l nsf l /fU)8 (15) 

where p* is the Coulomb sliding friction coefficient, p s is the 
maximum static friction coefficient, and 5 and are 
empirical parameters. 

To compute the electromagnetic engagement forces, f„\ and 
f n 2 , an electro-mechanical model of the EMA clutch is utilized 
(Eq. (16) and Fig. 5). When a command voltage, V c , is applied 
to an EMA clutch coil, magnetic flux is generated which 
travels axially through the clutch disk pack, across the 
nominal static air gap, h 0 , and the dynamic air gaps, h d , and 
then returns via the end plate. Since the coil circuit inductance 
is a function of the air gaps, the EMA coil current, /,., is 
computed from 


ApNyf2 ^ +R c I c =V c (t ) 
2 (h 0 +N d h d ) dt 


(16a) 


where N w is the number of coil windings, A p is the pole face 
area, p 0 is the free space permeability of air, and R c is the coil 
electrical resistance. Finally, the clutch engagement forces are 


fni = 


po A p N w 2 I c 2 
4(/? 0 +N d h d ) 2 


-k s h d \i= 1,2] 


(16b) 


where k s is the clutch disk separator spring stiffness. 
Furthermore, based on the magnetic saturation limit of the 
coils, the maximum EMA actuation force is with the clutch 
full engaged (i.e., with h d = 0) 


_ bo ApN w I c max~ r . _ , 

J m'max _ 7 ~TA V ~ 1 j2 J 


(17a) 


where T / sel is the dynamometer load torque setting and V d is a 
scaling voltage. Flere, T L , is essentially a constant resistive 
torque which approaches zero for small operating speeds. 

The clutch torques T c] and T c2 are computed based on the 
instantaneous condition of sticking or slipping as determined 
from the relative slip velocities between the clutch flanges and 
disks, = kl c jj - [;=1, 2], the axial engagement forces, /,' 

and f, 2 , and the clutch dry friction coefficient p s(fe . In the case 
of slipping, T ci , is computed as 

3 3 

Tci, slip = VstbfmNj MR °-, ~ Rd 2 ) sgn(n cfl - Q cdi ) (14) 

3 {R 0 ~ R d ) 


and the corresponding static breakaway clutch torque is 


Tkmax = 4{ \ ^ t* = 1 ^ ( 1 7b ) 

3 {R 0 --R d 2 ) 


With maximum current in the engaged condition 


lemax — fcmax — I sat ( 18 ) 

Rc 

where V c max is maximum applied EMA actuation voltage and 
I sat is the coil saturation limit. 
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Figure 5. — Electromagnetically actuated (EMA) clutch. 

To account for clutch stick-slip and clutch lock-up, four 
possible system states must be computed and checked at each 
model simulation time-step. Here, the four possible states are: 

I) Both clutches slipping (unconstrained) 

if | Q c f\ — Q-cdl I > £ I ^ cfl ~ ^ cd2 I ^ £ 
then T c] = T c]slip , T cl = T c i,slip (19a) 

Q/ — Q — §cf\ §cd2 §cf2 

II) Clutch 1 locked - Clutch 2 slipping 

if I Q c/1 - Q. cdl I < e & T c \ „ < T clmax 

then T cl =T clJI , T cl = T c2sUp (19b) 

q/7 = [<t>m ^ cd\ §cd2 §cf2 

III) Clutch 1 slipping - Clutch 2 locked 

if | £l c f2 ~ ^cd 2 I — £ & ^c2,777 — ^c2,max 

then T cl =T clslip , T c2 =T c2JII (19c) 

*\lll =[§m §cf\ <IWl §cd2 §lY 

IV) Fully locked (stalled) 



Figure 6. — Dual EMA Clutch/Gearbox system model. 


A schematic of the overall system analytical model based 
on Equations (8) to (19) is given in Figure 6. 

To compute the system response, the model is implemented 
using the Matlab/Simulink (The MathWorks, Inc.) environment 
where each sub-system (motor, EMA clutch coil, dual 
gearbox/shaft) are linked via their appropriate physical inputs 
and outputs. The model parameters used in the experimental 
validation tests are summarized in Tables II to IV. 


TABLE II.— DC MOTOR PARAMETERS 


Parameter Value 

Rotational inertia, J m 0.102 oz-in-sec 2 

Motor circuit resistance, R m 4.0 Ohms 

Armature inductance, L m 13.27 x le-3 H 

BackEMF constant, 1.1745 V/sec 

Torque constant, k mt 1.1745 Nm/A 

Proportional gain, k p 0.01 V-s 

Integral gain, kj 0.1 V 

Derivative gain, k d 0.001 V-s 2 

Set point speed, C2 sel 800 rpm 


else T cl - T cl IV , T c2 -T c2 IV 

(19d) 

q7F=[<i> m <lwi <\>cd 2 §lY TABLE III.— DYNAMOMETER LOAD PARAMETERS 


Parameter Value 

where £ is a small tolerance below which the slip speed is Rotational inertia, J L 96.06 oz-in-sec 2 

considered zero and the transmitted clutch torques T chII , T c2JII , Load torque set ting, T Uet 360 lb-in. 

Tcuv, and T cXIV are computed from a kinematic constrained Output shaft stiffness, k c2l , 11509N-m 

model with corresponding degree-of-freedom vector q//, q m 
and q lv for each locked condition. 
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TABLE IV.— EMA CLUTCH PARAMETERS 


Parameter 

Value 

Outer radius, R 0 

2.3125 in. 

Inner radius, R d 

0.5625 in. 

Number of disks, N d 

4 

Free space permeability, p 0 

4nxl0-7 N-A-2. 

Pole face area, A p 

4.43 in. 2 

Number of coil windings, N„ 

600 

Nominal Air gap, h 0 

0.005 in. 

EMA coil resistance, R, 

14 Ohms 

Maximum coil voltage, V cmax 

24 V 

Static friction coefficient, p, 

0.12 

Kinetic friction coefficient, jj. a . 

0.062* 

Stribeck velocity, 

0.001 

Stribeck exponent, 5 

1.0 

Max static torque, T ci ma , 

70 lb-l'f ' 

Max slipping torque, max(7,. /s | ip ) 

35 Ib-li" 


*@ 1800 rpra 
"with R= 24 V 


Chariot gearbox Torque and speed sensor Dynamometer 



Figure 7. — Chariot gearbox test setup. 


Apparatus 

Experimental validation tests were performed on the 
Chariot gearbox in the NASA Glenn gear noise test facility. A 
photograph of the test setup is shown in Figure 7. The Chariot 
gearbox was mounted on a custom built support table. The 
gearbox was attached to the support table through the same 
output mounting as used in the Chariot vehicle. Note that the 
gearbox was run horizontal in the test setup as compared to 
vertical in the Chariot vehicle. The output of the Chariot 
gearbox was connected to a commercially available 
torquemeter which measured gearbox output torque and speed. 
The output of the torquemeter was attached to the facility 
dynamometer which produced a resistance torque on the 
gearbox. The facility dynamometer was a water-cooled eddy 
current dynamic absorbing dynamometer with matching 
controller. 

A schematic of the instrumentation and control setup for the 
experiments is shown in Figure 8. The Chariot gearbox drive 
motor speed was controlled by a commercially available 
digital servo motor controller. Drive motor speed was 
manually adjusted by the operator to a preset value. Gearbox 
output torque was adjusted by the operator to a preset value 
using the facility dynamometer and controller. Low-speed and 
high-speed clutch engagement profiles were pre-programmed 
using a computer and an analog output card. Voltages from the 
analog output card for the low and high-speed clutches were 
routed to individual DC power amplifiers and then to the low- 
and high-speed clutches, respectively, thus allowing unlimited 
shifting potential. For each test, the transient Chariot gearbox 
output speed and torque, as well as the transient drive motor 
speed and current, were monitored using a facility data 
acquisition system. 


Chariot Gearbox 


Torquemeter 


Dynamometer 


Low 

speed 

clutch 


High 

speed 

clutch 


DC 

power 

amp 

— rz 


Chariot 

gearbox 

motor 

speed 

controller 


DC 

power 

amp 




Torque 


Analog 

output 

card 

1 


Speed / 
Current 


Speed 


Existing 
dynamometei 
controller 


Laptop 


Facility 

data 

acquisition 

system 


Monitor 


Laptop 


Figure 8. — Chariot gearbox instrumentation and control setup. 


Results and Discussion 

Two test results are presented. The first is a 1-sec low-to- 
high speed shifting profile (upshift). The second is a 
1-sec high-to-low speed shifting profile (downshift). For the 
upshift test, the low-speed clutch was engaged at the start of 
the test. The Chariot gearbox motor speed was set at 800 rpm, 
which produced a gearbox output speed of 50 rpm. The 
facility dynamometer was set at 400 lb-in., which produced a 
gearbox output torque of about 360 to 370 lb-in. as measured 
by the torquemeter. 

For the upshift case, the programmed 1-sec upshift profile 
was sent to the EMA clutches (Figs. 9(a) and (b)). Here, the 
low-speed clutch started fully engaged (24-V DC) and linearly 
decreased to fully disengaged (0-V DC) in 1 sec, while the 
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high-speed clutch started fully disengaged (0-V DC) and 
linearly increased to fully-engaged (24-V DC) in 1 sec. The 
start of the low-speed clutch disengagement occurred at the 
same time as the start of the high-speed clutch engagement. To 
simulate the upshift response, the same upshift command 
voltages (Figs. 9(a) and (b)) were sent to the EMA clutches in 
the analytical model during the simulation. 

The results of the upshift test are given in Figure 9. The 
start of the shifting occurred at 3 sec in the figure and 
completed at 4 sec. The Chariot gearbox output speed started 
to respond to the shift change at about 3.8 sec (Fig. 9(c)). A 
slight overshot occurred and the output speed settled to a 
steady value at about 5.5 sec. It should be noted that the 
overshot could have probably been avoided if further fine 
tuning of the motor controller would have been performed. 
However, that was not the objective of the test. The gearbox 
output torque exhibited a corresponding oscillation during the 
shifting process (Fig. 9(d)). This was probably due to the 
speed regime, where the output speed was at the low end of 
the dynamometer control range. The Chariot gearbox motor 
itself started to respond to the shift at around 3.4 sec as 
indicated by an increase in motor current (Fig. 9(g)). There 
was a slight drop of motor speed at the start of the output 
speed change, then a corresponding overshot as previously 
described (Fig. 9(h)). 

During test, the drive motor speed and drive motor current 
were measured from analog outputs of the motor digital servo 
controller. As a limitation of the controller software, speed and 
current could only be measured one at a time. Thus, the results 
from Figure 9 were actually from two separate inns. Identical 
initial conditions (motor speed and dynamometer torque) and 
identical shifting profiles were used for the two runs. When 
compared, the resultant transient Chariot gearbox output 
speeds (Fig. 9(c)) and torques (Fig. 9(d)) for these two runs 
were identical, thus validating the repeatability of the 

experiment. 

The corresponding simulation results for the upshift test are 
also given in Figure 9. The Chariot gearbox output speed 
simulation also predicted an overshot as with the experiments, 
with similar results in magnitude but slight phasing 

differences (Fig. 9(e)). The predicted gearbox output torque 
(Fig. 9(f)) and predicted drive motor current (Fig. 9(i)) also 
had oscillations during the shift similar to that from the 
experiments with similar magnitudes but slight phasing 
differences. The predicted drive motor speed did not have the 
initial drop in speed at the start of the shift as with the 
experiments, but it did have a similar oscillation correction 
after the shift (Fig. 9(j)). Overall, the predicted results 

matched quite well to the experiments, especially in 

magnitude of speed, torque, and current. It should be noted 
that the exact drive motor controller-speed control parameters 
were proprietary with the manufacturer. Thus, relative PID 
gain values could be derived from the control software but 
exact values used in the simulation had to be assumed. The 


control parameters of the drive motor had a major impact on 
the results, and thus could have lead to the discrepancies. 

For the downshift test, the high-speed clutch was engaged 
at the start of the test. The Chariot gearbox motor speed was 
set at 800 rpm, which produced a gearbox output speed of 
200 rpm. The facility dynamometer was set at 400 lb-in., 
which produced a measured gearbox output torque of about 
360 to 370 lb-in. The programmed 1-sec downshift profile was 
then sent to the EMA clutches (Figs. 10(a) and (b)). For this 
profile, the low-speed clutch started fully disengaged 
(0-V DC) and linearly increased to fully engaged (24-V DC) 
in 1 sec. The high-speed clutch started fully engaged 
(24-V DC) and linearly decreased to fully disengaged 
(0-V DC) in 1 sec. The start of disengagement for the high- 
speed clutch occurred at the same point in time as the start of 
engagement for the low-speed clutch. To simulate the 
downshift response, the downshift command voltages 
(Figs. 10(a) and (b)) were sent to the EMA clutches in the 
analytical model during the simulation. 

The results of the downshift test are given in Figure 10. As 
with the upshift tests, the start of the shifting occurred at 3 sec 
in the figure and completed at 4 sec. The Chariot gearbox 
output speed started to respond to the shift change at about 
3.3 sec (Fig. 10(c)). A very slight undershot occurred and the 
output speed quickly settled to a steady value at about 3.8 sec. 
The gearbox output torque had a significant drop during the 
shifting process, with the torque even changing directions 
(positive to negative in Fig. 10(d)). This was probably due to 
the dynamometer control during the speed decrease. The 
Chariot gearbox motor itself started to respond to the shift at 
about the same time as the gearbox output speed decrease 
(3.3 sec) as indicated by a slight increase in motor current 
(Fig. 10(g)). After that, there was a significant momentary 
decrease in motor current (down to zero) then an increase with 
an oscillation, then back to a steady-state value. There was a 
slight increase of motor speed at the start of the output speed 
change, then a corresponding overshot and return to a steady- 
state value (Fig. 10(h)). 

The corresponding simulation results for the downshift test 
are also given in Figure 10. The Chariot gearbox output speed 
simulation also predicted a rather quick drop in speed to the 
low-speed steady-state value with a very small oscillation 
(Fig. 10(e)). The predicted gearbox output torque (Fig. 10(f)) 
also had a significant drop during the shifting process and 
changed direction as in the experiments. The predicted torque 
exhibited oscillations during the torque drop. The predicted 
motor current also showed an increase, sharp decrease to zero, 
then increase to the steady-state value (Fig. 10(i)). The initial 
and final steady-state motor currents of the predictions were 
slightly lower than the measured results. Finally, the predicted 
drive motor speed had a very similar oscillation during shift as 
that for the experiments (Fig. 10(j)). Again, the predicted 
results matched quite well to the experiments, in magnitude 
and even phase, for speed, torque, and current. 
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Figure 9. — Comparison of predicted and measured transients for a 1-sec upshift profile. 
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Conclusions 

Dynamic simulations and experimental validation tests 
were performed on a two-stage, two-speed gearbox as part of 
the drive system research activities of the NASA Fundamental 
Aeronautics Subsonics Rotary Wing Project. A dynamic 
model of the system was created which included a DC electric 
motor with proportional-integral-derivative (PID) speed 
control, a two-speed gearbox with dual electro-magnetically- 
actuated clutches, and an eddy current dynamometer. A six 
degree-of-freedom model of the gearbox accounted for the 
system torsional dynamics and included gear, clutch, shaft, 
and load inertias as well as shaft flexibilities and a dry clutch 
stick-slip friction model. The simulation was used to predict 
gearbox performance during shifting scenarios. Experimental 
validation tests were performed on the gearbox in the NASA 
Glenn gear noise test facility. Gearbox output speed and 
torque as well as drive motor speed and current were 
compared to those from the analytical predictions. The 
following conclusions were obtained: 

1) Overall, the predicted transient results matched quite 
well to the experiments, especially in magnitude of speed, 
torque, and current. 

2) The control parameters of the drive motor had a major 
impact on the results, and thus, could have lead to the 
discrepancies since exact PID values used in the simulation 
had to be assumed. 

3) The transient Chariot gearbox output speeds and torques 
from the experiments were extremely repeatability. 
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